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We investigate the persistent current of a ring with an in-line quantum dot capacitively coupled 
to an external circuit. Of special interest is the magnitude of the persistent current as a function 
of the external impedance in the zero temperature limit when the only fluctuations in the external 
circuit are zero-point fluctuations. These are time-dependent fluctuations which polarize the ring- 
dot structure and we discuss in detail the contribution of displacement currents to the persistent 
current. We have earlier discussed an exact solution for the persistent current and its fluctuations 
based on a Bethe ansatz. In this work, we emphasize a physically more intuitive approach using 
a Langevin description of the external circuit. This approach is limited to weak coupling between 
the ring and the external circuit. We show that the zero temperature persistent current obtained 
in this approach is consistent with the persistent current calculated from a Bethe ansatz solution. 
In the absence of coupling our system is a two level system consisting of the ground state and the 
first excited state. In the presence of coupling we investigate the projection of the actual state on 
the ground state and the first exited state of the decoupled ring. With each of these projections we 
can associate a phase diffusion time. In the zero temperature limit we find that the phase diffusion 
time of the excited state projection saturates, whereas the phase diffusion time of the ground state 
projection diverges. 

PACS numbers: 73.23.Ra, 73.23.Hk, 71.27.+a 



I. INTRODUCTION 

An interesting aspect of quantum systems is the fact 
that even at zero temperature there are fluctuations 
which manifest themselves in observables which do not 
commute with the Hamiltonian of the system. In meso- 
scopic systems phase coherence plays an essential role, 
and it is thus important to ask to what extent coherence 
is affected by zero-point fluctuations. To investigate this 
question we consider the ground state of a normal meso- 
scopic ring threaded by an Aharonov-Bohm flux and ca- 
pacitively coupled to an external circuit. Such a ring ex- 
hibits a persistent current which is a direct consequence 
of phase coherent electron motion over distances large 
compared to the ring circumference. At zero tempera- 
ture such a ring interacts with an external circuit (see 
Fig. |IJ) only due to zero-point fluctuations. More pre- 
cisely, the external circuit can, through the generation of 
zero-point voltage fluctuations, induce polarization fluc- 
tuations in the ring which in turn affect the magnitude 
of the persistent current. The source of the voltage fluc- 
tuations in the external circuit are the resistive elements. 
We can thus ask: How does the persistent current of the 
ring depend on the resistive properties of the external cir- 
cuit? For the system shown in Fig. [j], Ref. [j] provided an 
answer by mapping a simple model of a ring with a quan- 
tum dot and external circuit on the anisotropic Kondo 
model and using the known Bethe ansatz solutions of 
this problem.0 The purpose of this work is to consider the 
same model and to provide a discussion which is phys- 
ically more transparent. The discussion given below is, 
however, limited to the case of weak coupling between 



the ring and the external circuit. 




FIG. I. Ring with an in-line dot subject to a flux "I> and 
capacitively coupled to an external impedance Z ex t- 

Investigations of the coherence properties of the 
ground state of mesoscopic structures are rare and con- 
cern mainly superconductors. Hekking and GlazmarfcJ in- 
vestigate a thin superconducting loop in an.jelectromag- 
netic environment; Oshikawa and ZagoskinQ consider a 
superconducting grain coupled to a normal wire. In con- 
trast, investigations on dephasing in systems driven out 
of the ground state are numerous and there exists a con- 
siderable literature. In particular, zero-point fluctuations 
have been of interest, following recent experimental work 
by Mohanty, Jariwala, and WebbO on weak localization 
in metallic diffusive wires. These experiments probe a 
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transport state and not the ground state of the system. 
In the experiments of Mohanty and Webb it is the weak- 
localization effect in the linear response conductance of 
the system which is of interest. Subsequently to these 
experiments the role of zero-point fluctuations, in weak 
localization effects has been very much debatedEITJ and a 
number of works have appeared which suggest that the 
experiments can possibly be better explained without in- 
voking zero-point fluctuations.!!^ Weak-localization is a 
quantum effect which survives ensemble averaging and 
thus the dephasing rates which count are specific to the 
fact that we deal with time-reversed trajectories and that 
an ensemble average has to be performed. The absence 
of an effect of zero-point fluctuations found in Refs. |t],|| 
might give rise to the mistaken impression that zero-point 
fluctuations are quite generally without any effect on the 
coherence properties of a system. The example presented 
here shows that zero-point fluctuations are clearly impor- 
tant for the coherence properties of the system, although 
perhaps not for weak localization phenomena. Below we 
consider a specific ring (not an ensemble) and ask how 
the maximum amplitude of the persistent current is af- 
fected by zero-point fluctuations. 

The effect of a thermal bath on the persistent, cur- 
rent hasJaeeia discussed by Landauer and BiittikerU and 
Biittikeni3'E3 within a Debye relaxation approach. In 
this approach, however, the bath affects essentially only 
the population of different states, but not the electronic 
states of the ring itself. As a consequence, at zero tem- 
perature, this model exhibits no effect due to zero point 
fluctuations: the magnitude of the persistent current is 
independent of the coupling strength to the bath. Still 
a different model, introduced by Buttiker,E-3 considers a 
ring coupled via a side branch to a normal electron reser- 
voir. Due to the connection between reservoir and ring 
a carrier in the ring eventually escapes to the reservoir 
and is replaced with an incident carrier with a phase 
which is unrelated to that of the escaping carrier. This 
model predicts even at zero temperature an amplitude 
of the pexsastent current which depends on the coupling 
strength.L3~Ej This effect is not due to zero-point fluctu- 
ations but results from the exchange of carriers between 
the reservoir and the ring. If the ring is coupled to a 
side branch of finite length, the side branch can nev- 
ertheless generate effects which are similar to that of a 
reservoir, especially if only ensemble averaged quantities 
are considered. This is correct only if the side branch has 
a charging energy which is weak compared to the level 
spacing.113 If the charging energy is large compared to 
the level spacing the side branch hasjio effect on the en- 
semble averaged persistent current &B It is thus interest- 
ing to ask, whether there exist models which are strictly 
canonical (without carrier exchange with a reservoir) and 
for which the sample specific persistent current depends 
nevertheless on the properties of the bath. The model in- 
vestigated in this work examines a ring which is coupled 
to the bath only via the long range Coulomb force. 

The mesoscopic ring which we consider (see Fig. [l]) 



is divided into two regions by tunnelijig^barriers. It is 
a ring with an in-line quantum dot.E9~Ej. This model 
allows a simple characterization of the electrostatic po- 
tential in terms of only two variables Ud (for the dot) 
and U a (for the arm), and by the charges Qd and Q a . 
The two regions, the dot and the arm of the ring, are 
coupled via capacitors C\ and C2 to the external circuit. 
The external circuit is described by its impedance Z ext . 
The potential at the capacitor C\ is denoted by Vo, the 
charge by Qo- Likewise, we write the potential and the 
charge on the capacitor C2 as Voo and Qoo- 



II. LANGEVIN EQUATION 

To be specific, we consider a purely resistive external 
impedance, Z ext — R. The fluctuations generated by this 
resistor can be represented as a noise source in parallel 
with the resistor as shown in Fig. ||. The current 5V/R 
through the resistor must be equal to the sum of the cur- 
rent of the noise source 6I(t) and of the current through 
the ring-dot structure. In terms of the potential differ- 
ence SV = Vq — Voo across the ring-dot structure we find 
that the coupling between the ring-dot structure and the 
external circuit is described by 



SV/R = -C SV 



Co 

Ci 



Q d + SI(t). 



(1) 



Here we introduce the capacitance Cq 1 = C~[ + C~ x 
which is the series capacitance of the internal capac- 
itance Ci = Cl + Cr and the external capacitance 
C~ 1 = C^ 1 + C% , and the charge on the dot Qd- A 
detailed consideration of the circuit equations leading to 
Eq. (0) is given in Appendix |X| The noise source in par- 
allel to the resistor generates a current noise spectrum, 
(SI{uj)5I(uj')) = 2ttS{w + w')Sh[u), given by 



„ . . huj , / Tvjj \ 

S n {uj) = — coth 

17 R \2kT 



(2) 



where k is the Boltzmann constant, and T is the tem- 
perature. To complete the description of this system, 
we need to investigate the dynamics of the ring in the 
presence of the fluctuating external voltage 5V(t). 
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FIG. 2. The external impedance Z ext (uj) has been replaced 
by a resistor R with a noise source in parallel in order to take 
into account quantum fluctuations. 

First, we consider the effect of the external circuit in 
linear response. Since we are only interested in the weak 
coupling limit, we can choose the ground state of the 
ring-dot system in the absence of an external circuit as a 
state around which we can expand. What we need is the 
linear response relation between 5Qd{t) (the deviation of 
the charge away from a reference state) as a function of 
the external voltage SV(t). Note that both of these vari- 
ables are in fact operators. Here, we proceed as in the 
random phase approximation and use a linear response 
function to describe the connection between these two 
quantities. In linear response we can Fourier transform 
SQd(t) and SV(t). The linear response coefficient gives 
the increment of the charge on the dot in response to a 
variation of an external voltage: It is thus a dynamic ca- 
pacitance which we denote by Cd{oj) — 6Qd{w)/8V(u)). 
Similarly, the relation between the charge SQq piled up 
on the external capacitor C\ and the applied voltage is 
determined by the overall capacitance of the ring dot- 
structure vis-a-vis the external circuit, and is denoted 
by Cp(w) = SQ (uj)/SV(lu). From circuit theory (see 
Appendix [A]) it follows that these two capacitances are 
related 



Co - ^QH 



(3) 



The noise source (see Fig. ||) sees a total impedance Z{uj) 
which consists of this capacitance and the resistance 



Z{oj) 



1 

R 



(4) 



With these response coefficients Eq. ([j]) leads to 

8V{uj) = Z{u))8I{u) = [1 - iuRC^w)}' 1 R5I{u) (5) 
and the spectrum of the external voltage is thus 



S vv (uj) 



Rhu coth(^) 
[l+a,2(i?C M M)2] 



(6) 



Our next task is now to find explicit expressions for the 
capacitances C^(lo) and Cd{ui). 



A. The dynamics of the ring 

The dynamics of the charge on the dot, determined by 
Qd, is in the Heisenberg picture given by 



d - h ~ 

TTiQd — ~[Rrinq , Q d] , 

at I 



(7) 



where H r i ng is the Hamiltonian of the ring in an applied 
external potential 6V. 



We assume that the tunneling amplitudes are much 
smaller in magnitude than the level spacing in the dot 
and the level spacing in the arm. Moreover, we assume 
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the electrons to be spinless. As discussed in Rcfs. |o 
spin singlet appears in the case of electrons with spin, and 
the tunneling amplitude is enhanced by a factor of v^i 
Our choice of parameters excludes the Kondo effect,Ej 
also if the electron spin is taken into account. As a mat- 
ter of fact, the Kondo effect appears only at a tunnel- 
ing amplitude-.comparable to or larger than the mean 
level spacingrJ We note that the number of charge car- 
riers in |thc_ring is conserved. Following Buttikcr and 
StaffordpJ'EJ we consider hybridization between the top- 
most occupied electron level in the arm and the lowest 
unoccupied electron level in the dot, e a M and e d [N+i) 
only. To simplify the notation, we denote the tunnel- 
ing amplitudes between the levels e a M and e^jv+i) by 
th for tunneling through the left junction and by tn for 
tunneling through the right one, and introduce the total 
tunneling amplitude 



t\ ± 2thtR cos 



2tt$ 



(8) 



where $0 = hc/e is the flux quantum, and the sign in 
front of the cosine depends on the parity of the num- 
ber of electrons in the ring. The sign is positive, if the 
number of electrons is odd, and negative if the number is 
even. In the model introduced above, we are considering 
only two states, namely the topmost electron being in 
the dot, which we represent by the vector (1, 0), and the 
topmost electron being in the arm of the ring, written as 
(0, 1). The dynamics of these two states is described by 
a time dependent Hamiltonian 



H 



he(t) 



hA 



hv{t) . 



ring 



(9) 



where a z and a x are Pauli matrices, and 1 is the identity 
matrix. The prefactor Tiv(t) = —Co 5V 2 (t) is a global 
shift in energy. We split up the detuning e(t) into a time 
independent and a time dependent part e(t) = So + Ss(t), 
with 



e 



e 2 (iV-iV + + l/2) C„ r 

— --c Vo 



c 



c 



{Qd* — Qdo) 



(10) 

(11) 



The effective background charge on the dot is Qdo = 
eN + + C e Vo, where the first term is a built in background 
charge and the second term can be externally controlled 
by applying a static voltage across the ring-dot struc- 
ture. If coherent tunneling is neglected the ground states 
of the ring with N and N + 1 carriers are degenerate if 
the polarization charge is equal to 
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Qd* = e \N + M + j (e d{N+1) - e aM ) ■ (12) 

Thus Qdo = Qd* is the (classical) condition at which the 
Coulomb blockade is lifted. Quantum mechanically the 
state in the ring and the state in the arm of the ring are in 
resonance when Qdo — Qd* an d the persistent current ex- 
hibits a peak. Our simple two-level picture is applicable 
when £o is small, that is, in the vicinity of a resonance. 

We want to find the time evolution of a state ip(t), 
which we write as 



/cos § e^/ 2 
Vsin | e -^/ 2 



(13) 



with 9, ip and \ real. This is the most general form 
of a normalized complex vector in two dimensions. In 
terms of 9, <p and the global phase x> the time dependent 
Schrodinger equation reads 



<P 



X 



— e — Se(t) — A cot 9 cos ip, 
-A sin^, 

COS If 



Ao 



sin 9 



(14) 
(15) 
(16) 



Note that Eq. ( |l6|) describing the dynamics of the global 
phase x is special in the sense that it is driven by <p and 
9 but it has no back-effect onto them. Moreover, \ is 
irrelevant for expectation values, like the persistent cur- 
rent or the charge on the dot. It plays an important role, 
however, in the discussion of phase diffusion times, see 
Sec. [TV|. The system of the two equations for <p and 9, 
Eqs.^Tj, [L5|), is closed by Eq. (|l|), 



Qd=j t m)\Qd\m), 



with ip(t) defined above, and 



1 




■eN. 



+ • 



Wc find 



Qd = e-m)\ 



W)) = 



e 

--smt 



Using this result and Eq. ([[J), gives 



SV = -CqRSV 



Coe 
Ci 2 



sin9 9 + RSI(t). 



(17) 



(18) 



(19) 



(20) 



Eqs. ( |l~4] , |l5| ) and ( p(i| ) form a closed system of equations 
in which the external circuit is incorporated in terms of a 
fluctuating current SI(t) and of an ohmic resistor R. In 
the next section, we investigate Eqs. ( |l4| , [l5| ) and (|2(]) to 
find the effect of zero-point fluctuations on the persistent 
current of the ring. 



B. Expansion around a stationary state 

First, let us discuss the stationary states of the system 
of differential equations, Eqs. ([H], |l5|) and ( |20| ) in the 
absence of the noise term 51 (t). We take < <p < 2ir 
and < 9 < ir. This gives immediately sin ip = and 
consequently a stationary state has ip = tpo, with tpo = 
or ipo = it. With this it is easy to show that in the 
stationary state we must have 9 = 9o, with 



cot ( 



±- 



"A ' 



(21) 



The lower sign applies for (po = 0. This is the ground 
state for the ring-dot system at fixed e(t) = eg, and the 
upper sign holds for po = tt. The energy of the ground 
state is — HQ /2, thus the global phase is xo{t) = ^M- 
Here 



n 2 

5 2q 



(22) 



is the resonance frequency of the (decoupled) two-level 
system. We also introduce the "classical" relaxation time 
trc = RCq, and a relaxation rate 



A 2 



(23) 



which as we shall see in Sec. [IE is a relaxation rate due 



to the coupling of the ring-dot system to the external 
circuit. Here a is a dimensionless coupling constant 

(24) 



R 
Rk 



Co 



with Rk = h/e 2 the quantum of resistance. We will 
immediately see the usefulness of these definitions. 

Now, we switch on the noise 5I(t). We seek <p(t), 6(t), 
x(t) and SV(t) in linear order in the noise current SI (t). 
We expand <p(t) and 9(t) to first order around the ground 
state, <p = and 9 = Oq. For Sp(t) = p(t) — tpo, 
66(t) = 9{t) — 9q, we find in Fourier space, 



Q2 

iivSp = -Se + -T-S6, 
A 



iu;S9 



jSe = 



-A S(p, 
1 



trc 



-Se - r<5^ 



I4« 



(25) 
(26) 
(27) 



We also expand the global phase x(f) around its evo- 
lution in the ground state Xo(t) = and define 
^xW = x(^) — Xo(t)- I n Fourier space, Eq. ( |l6| ) becomes 



iivSx = ^n-r— t 
A 



(28) 



We note that there is no effect of the global shift in en- 
ergy, hv(t), as it is quadratic in the voltage SV, and we 
are only interested in effects up to linear order in SV. 

In the following section, we evaluate the linear response 
of the ring, described by Eqs. (^5|, ||^) to an applied ex- 
ternal potential 5V{oj), giving the frequency dependent 
capacitance C^(lo). 
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C. Capacitance of the ring and impedance 

We evaluate 8Qd{t-o) to first order in SV (u), using 
Eqs. (|25|, ||^), which gives the dynamic capacitance 
C d {u) = 5Q d (u)/5V(u), 

e 2 /(2C % ) 



C d (to) = C - 



A 2 



n 2 ' 



(29) 



The frequency dependent capacitance as seen from the 
external circuit reads, see Eq. (j^), 



C„(w) =Q)1 



e 2 /(2C.) 



A 2 



fi 2 
"0 



(30) 



The frequency dependent capacitance C^ contains in ad- 
dition to the geometrical capacitance Cq a term which 
arises from the dynamic polarizability of the ring. In 
the low frequency regime the polarizability enhances the 
capacitance, 



Cm(o) = 0)1 + 



e 2 /(2C») A 2 



n 2 



Co 



2e 2 \t ± f 



{{e(Q d o-Qd,)/C} 2 + 4|t±| 2 } 



3/2 ' 



(31) 



The zero frequency capacitance has already, been ob- 
tained by Stafford and one of the authors-E^IO The zero- 
frequency (electro-chemical) capacitance Eq. ( |3l| ) can be 
obtained as a-second order derivative of the grand canon- 
ical potential 2 .!! with respect to the voltage applied across 
the ring dot-system or as a dynamic response to a slowly 
varying voltage across the ring-dot system. Since the tun- 
neling amplitude Aq depends periodically on the mag- 



netic flux $, see Eq. (| 



Cn(u)) can be modulated by 



varying The electro-chemical capacitance exhibits a 
peak as the bias Vq is varied. The peak shows up at reso- 
nance, Qdo — Qd*, corresponding to So = 0. At the same 
point, there is also a peak in the persistent current. A re- 
cent experiment by Deblock et alxHi demonstrated indeed 
a flux dependent polarizability of mesoscopic rings. 

At high frequencies the polarization cannot follow the 
external voltage and the capacitance Eq. (|3^) is smaller 
than the geometrical capacitance. At u = flo the capac- 
itance exhibits a resonance characteristic for dielectric 
functions. The ring by itself has no damping mechanism, 
therefore the resonance shows up as a pole. The damping 
is provided by the dissipation in the external circuit. 

From Eqs. (0) and (Q, we find 



1 



lw 2 -!lg - iuT R c (uj 2 — 0§ — T/t rc ) 



Z(uj) R 



O 2 



(32) 



This expression for the impedance seen by the noise 
source contains all the information needed to calculate 
the various spectral densities. 



89(uj) and their cross-correlations can all be expressed 
in terms of the current noise spectral density. We have 
See{u) = (AI/uj 2 )S vv (uj), S v e{oj) = -i(Ao/uj)S vv (uj), 
and 



S v Jw) 



Co 

a 



" 2 \Z(co)\ 



ij {uj 2 ~%) 



(34) 



We will relate the reduction of the persistent current to 
uj) in Sec. III. While the spectral densities for Sip 
at zero frequency, the spectral den- 



and 



are finite 
sity of the global phase Sx is not. 



We have S xx (uj) 



sec 



Eq. (t 



-pole in the vicinity of uj 



|), therefore it exhibits a 
for finite temperatures 



S xx (uj) ~ 2na 



2kT el 1 



(uj - 0), 



(35) 



whereas in the quantum limit T 
is reduced by one, 



: the order of the pole 



D. Spectral densities 



We solve Eqs. (|25|)-(|27|) for Sip by eliminating 89 and 
Ss, and obtain after some algebra 



S(p = e C ° ^~ iuJ "> Z ^ §j 



n 2 



(33) 



It follows immediately that (Sip(u>)) = (56(u>)) = and 
that (Se(uj)) = 0. The spectral densities of Sip(ui), 



S xx (uj) ~ 2na 



- 2 1 



(36) 



see Fig. [|. These two last results do not affect the per- 
sistent current, but they are of great importance in the 
discussion of the rates of phase diffusion in Sec. IV. 
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FIG. 3. The spectral function S xx (ui) at finite temperature 
kT — lOTifio (solid line) and at zero temperature (dashed 
line). The left y-axis carries the scale for the finite temper- 
ature curve, the right one carries the scale for the zero tem- 
perature curve. The parameters are the same as in Fig. ^, in 
addition eo = fio/v^- 



E. Weak coupling 

Within the approach out-lined above, we can expect to 
reproduce the exact result for the reduction of the persis- 
tent current only in the limit of weak coupling between 
the ring-dot system and the external circuit. This means 
that the potential fluctuations of the system and the ex- 



ternal circuit are only weakly coupled. Thus weak cou- 
pling between the ring and the external circuit is achieved 
by letting the capacitances coupling the two systems be- 
come very small,Ea C\, C2 — ^ 0, for a fixed resistance R. 
This entails Co — ► 0, thus the coupling str eng th a be- 
tween the ring and the external circuit, Eq. (|24|), is small 
against 1. Note that in this approximation the RC-time 
trc becomes very small as well. 

In the following paragraph, we calculate the poles of 
Z(uj), in the approximation Cq — > 0. The small parame- 
ters in this case are trc oc Co and T cx Cq . The equation 
= has one solution behaving as Cq . To lead- 
ing order in Co, it reads 



0J3 



1 

trc 



(37) 



Thus in the weak coupling limit charge relaxation across 
the mesoscopic system becomes instantaneous. The other 
two solutions, u>±, are of zeroth order in Co- Up to and 
including corrections of order Cq , they are 



uj ± = ±n - i~, 



(38) 



with a relaxation rate T introduced in Eq. ( |23| ) . The rate 
r describes the relaxation of perturbations with a fre- 
quency close to the eigenfrequency of the isolated ring- 
dot system. In Sec. IV it is shown that F is also related 
to a phase diffusion time. The denominator in Eq. (|3~i| ) 
can now be written as 



\Z(w)\* 



1 RC 

R 2 



' RC 



(w - n ) 2 + (r/2) 2 ( w + n ) 2 + (r/2) 2 



For frequencies, ujtrc ^ 1> the denominator takes the form 



or 



\Z(coW 



1 

fi2 



(w - n ) 2 + (r/2) 2 ( w + n Q ) 2 + (r/2) 2 



(39) 



(40) 



Since in the quantum limit it is always necessary to introduce a cut-off, Eq. (|40|) is valid over the entire range of 
frequencies of interest, if only 1/trc is larger than the cut-off frequency. The spectral density Eq. ([34]), reads 

therefore 



2irauj 3 coth ^ 



(w - r> ) 2 + (r/2) 2 (lu + n ) 2 + (r/2) 2 



(41) 



where the coupling strength a between the ring and the 
external circuit is defined in Eq. (gj). Note that S vv (uj) 
goes to zero at the origin uj = like \lu\ 3 in the extreme 
quantum case (T = 0), and like to 1 for finite tempera- 
tures, see Fig. m For |ct>| ^> fig, it behaves like for 
small temperatures and like lo~ 2 for large temperatures. 
Analogous formulae hold for Sgg(ui) and for S v g(uj). We 
mention that the spectral density Seg(uj) behaves for 



small u) like \ui\ at T = and goes to a constant for 
finite T. All the three spectral densities exhibit peaks of 
width r at ±Oq- 
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The spectral density of the global phase S xx (to) , Fig. ||, 
exhibits peaks at to — ±f2 as well. More important, how- 
ever, is the fact that it has a pole at to = 0, as already 
pointed out in the previous section. 



III. PERSISTENT CURRENT 



In the following, we calculate the persistent current as 
a quantum and statistical average of the operator of the 
circulating current 



FIG. 4. The spectral function S^ v (to) at finite temperature 
kT — 107if2o (solid line) and at zero temperature (dashed 
line). The left y-axis carries the scale for the finite temper- 
ature curve, the right one carries the scale for the zero tem- 
perature curve. The parameters are R = 5Rk, Cq = O.ICi, 
A = 0.5f2 and T = 0.2fi o . 



where J is given by 



J 
J* 



(42) 



J 



t\ +t 2 R ± 2t L t R cos 



2vr$ 



-1/2 



±t L t R sin ■ 



2?r$ 



Cb 9 CY 9 C R — Ct, 27r$ 
— t ? -tl ± — -t L t R cos 

Ci d R a L R $n 



(43) 



Here the first term is a pure particle current contribution, 
whereas the second term is a consequence of interactions. 
For the average persistent current of interest here, it is 
only the first term which contributes. A derivation of 
Eq. (|^) and a discussion of the relationship between I c 
and the magnetization is given in Appendix |b| 

The expectation value of the persistent current for the 
state given in Eq. ( |l3| ) reads 

7(f) = (^(t)|J c |V(t)) - ^RBjsmBe-*. (44) 

We are, however, interested in the statistically averaged 
of the persistent current (/(£)), given by 

(W«)|/<#(*)» - ^Re (J{sm8e-*)) , (45) 

where the double bracket indicates a quantum and sta- 
tistical average. Therefore, we have to calculate the cor- 
relator (sm9 e~ t(p ) . First, we observe that there are no 
correlations between Sip and S9, 



(8<p(t)86{t)) = 0, 



(46) 



as the spectral density S v $(to) is an odd function of to. 
Thus we have 



-Vp\ 



-vp\ 



(47) 



SI. Second, it can be shown that the second order cor- 
rections to (p and 9 vanish on the average. Finally, we 
assume that the correlations in Sip are Gaussian, which 
allows us to write 



l e - i( P(t)\ _ e -»VD/ e -***»(t)\ 



exp 



. (48) 



where we have used that ^0 = 0. In the weak coupling 
limit, see Sec. [II E| , and in the extreme quantum limit, 
T = 0, we find for the time averaged mean square fluc- 
tuations 



(VW) = / —S vv (co) 



2a In— , 



(49) 



where the cut-off frequency to c is taken to be larger 
than all frequency scales (except t rc ), and a is the 
coupling strength between ring and external circuit, see 
Eq. J24|). In the limit to c > Q , we can neglect (S9 2 (t)) — 
J^ a (duj/TT)See(io) against (S(p 2 (t)). We insert (5tp 2 (t)} 
and sin 0o = Ao/Slo into (sin^e -1 ^), and observe that 



2 J 2 a$ ■ 



(50) 



Put together, we obtain the noise averaged current in the 
ring 



to second order in 51, that is, the probability distribu- 
tions for ip and 9 are decoupled up to second order in 



(Ht)) 



he 8A A / n 



(51) 
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Eq. ( plj ) is a key result of this work. For a <C 1, cor- 
responding to weak coupling between the ring and the 
external circuit, the power law for the persistent current 
obtained in Eq. (|5l"l), as well as the exponent a, Eq. (24), 
are the same as the one obtained when the external cir- 
cuit i&,treated quantum mechanically as well. In a recent 
work,EJ the authors in collaboration with Ponomarenko 
have shown that if the external circuit is represented by 
a transmission line, the persistent current at resonance, 
e = 0, has for a < 1 the power law behavior 



7(e = 0) oc I — 

UJ r 



(52) 



For a very small coupling parameter, aCl, the Bethe 
ansatz result Eq. ((52|), goes over to the power law of 
Eq. (jjill). Thus the simplified discussion presented here 
leads at least in the weak coupling limit to the same re- 
sult as the one obtained in Ref. [jj 

We emphasize that the persistent current is a property 
of the ground state of a system. In our case, the per- 
sistent current is, however, carried by only a part of the 
system. Due to the coupling to the external circuit this 
subsystem is subject to fluctuations which even at zero 
temperatures suppress the persistent current. If we keep 
the capacitances fixed, Eqs. ( |5l| ) and ( f^ ) give a per- 
sistent current which decreases with increasing external 
resistance R. The spectral density of the voltage fluctu- 
ations across the ring dot system increases with R (see 
Eq- (0))- We next characterize the fluctuations of the 
ring-dot subsystem in more detail. 



IV. PHASE DIFFUSION TIMES 

In this work we are concerned with the effect of an ex- 
ternal circuit on the ground state properties of the system 
and in particular on the persistent current of the ring- 
dot subsystem. In transport experiments which reveal 
phase-coherent contributions or in experiments which in- 
vestigate the evolution of an initial state, we can describe 
the effect of an environment with the help of dephasing 
rates. On the other hand, it is not a priori clear that the 
coherence properties of the ground state of a system, can 
be equally expressed in terms of dephasing rates. After 
all the persistent current investigated here exists quite in- 
dependently on how long we let the system evolve. Here 
we consider the state in the ring and investigate its evo- 
lution away from an initial state. Below we show that 
this evolution is diffusive at least for times short enough 
such that we remain within the limit of validity of the 
calculation. 



We have seen in Sec. II E that the dynamics of the ring 
strongly affects the spectral densities 5 w (o;) and Sgg(u>) 
in the vicinity of the characteristic frequency f^o- The 
characteristic frequency describes the free dynamics of a 
ring which is disconnected from the external circuit, de- 
scribed by the Hamiltonian H r i ng , Eq. (^), without the 
time dependent terms 5e(t) and v(t), 



H = — a z 



(53) 



This Hamiltonian can be viewed as describing a spin in a 
magnetic field of strength proportional to (eq + A 2 ,) 1 / 2 = 
£!o, forming an angle 9q wit h th e z-axis. The angle 
6>o has been defined in Sec. |ll B| and is related to eq 
and Ao by cot#o = — eo/Ao- The eigenstates of the 
time independent Hamiltonian Hq are the ground state 
ip_ = (cos6*o/2, sin6*o/2), with eigenvalue —TiQ,q/2, and 
the excited state = {— sin^o/2, cos#o/2) with eigen- 
value hfly/2. The evolution due to H of a spin prepared 
in a state which is not an eigenstate corresponds to a 
rotation of the expectation value of the spin about the 
direction of the magnetic field with a frequency Qq. Let 
us now return to the full problem in which polarization 
fluctuations modify the free evolution of the decoupled 
system. In the weak coupling limit, we are interested 
in the time evolution which is long compared to Qq 1 - 
Therefore, we switch to the "rotating frame" picture, 
where the wave function ip(t) of the ring, Eq. ([l3|), reads 
ipn(t) = exp(iH t/h)ip(t). We consider the projection 
c±(t) = (i'±\tpR(t)) of the wave function ipn(t) onto the 
states ip±. If the state ip(t) evolves under the influence 
of the time independent Hamiltonian Hq, the projections 
c±(t) are constant in time. The moduli of the projections, 
averaged over the noise, (\c±(t)\ 2 ) are also independent 
of time if the evolution is determined by H r i ng which 
includes the fluctuating potential. The phase of the pro- 
jections c±(t), however, shows diffusive behavior. On the 
average, we have (\c±(t) — c±(0)| 2 ) ~ t/r± for sufficiently 
long times t. We shall see below that t_ is related to the 
phase diffusion of the global phase whereas t + is re- 
lated to the diffusion of the internal phases tp and 9, see 
Eq. (|l3|) . It is a known feature of two-level-systems that 
they exhibit two distinct dephasing times.Ea The differ- 
ence between r_ and r+ is particularly pronounced in 
the low temperature limit. The phase breaking time r_ 
diverges at zero temperature, whereas the time t+ satu- 
rates to a finite value for temperatures kT below MIq. 

We expand the wave function ipR to first order in 8e, 
that is, to first order in Sip, 89 and Sx, 



£ Slf(t) \ 

n 2 J 



89{t) .A d(p(t) 

~ l no~ 2 



(54) 



From this expression, we immediately obtain the projections c_(t) = 1 + i[Sx(t)/2 — (eo/flo)5ip(t)/2] and c + (t) = 
l/2(S9(t) - i(A /n )S<p(t)) exp(in t). As the averages (5ip 2 (t)), (S9 2 (t)) and (Sx 2 (t)) as well as all the "crossed" 
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averages of the type (8<p(t)86(t)) are constant in t, the moduli (\c±(t)\ 2 ) are indeed independent of time. The mean 
squared displacement of the projections c±(t) read respectively 



|c_(t)-c_(0)|- 



C+ (*)-c + (0)|- 



efcj 2 Ujt 

tt - sin — 

2vr 2 

2?T 2 



"0 

A 2 

+ fi ) + Yf2 S w( 



(55) 
(56) 



The long time behavior of Eq. ( |55| ) is dominated by the 
frequencies near u) = 0. The spectral density S vip (uj) van- 
ishes like oj 2 for finite temperatures or even like |w| 3 in 
the zero temperature limit. The spectral density S xx (lo), 
however has a pole at lu — 0. For finite temperatures the 
pole is of order io~ 2 , see Eq. (|35|), entailing a long time 
behavior of the type (|c_(i) — c-(0)| 2 ) ~ £/ t -j with a 
dephasing time 



n 2 



2?ra fcT eg 



(57) 



We stress again that r_ determines the diffusion time of 
the global phase x as it is given in terms of the spectral 
function S xx (u>). In the quantum limit T = 0, the phase 
diffusion time t_ diverges. Furthermore, it follows from 
Eq. ( p7[ ) that r_ is tunable, since £o may be varied by an 
external DC-bias. In particular, at resonance £o = 0, the 
phase diffusion time r_ diverges for any temperature. 

The long time behavior of Eq. ([56]) , on the other hand, 
is determined by the frequencies near Qq. In the vicinity 
of this characteristic frequency, S vv (tto + lu) as well as 
See(Qo + lu) show a lu~ 2 behavior at finite as well as at 
zero temperature, which is cut off by the relaxation rate T, 
defined in Eq. (|23| ) at very small frequencies lu ~ T <C • 
In summary, we have, for \ui\ -C fio 



2na Q coth |^ 

' w)ra w 2 + (r/2)2 



(58) 



and ^^(rio 



(Ag/^)5 w (fi + w). The time evo- 



lution of Eq. ( p6|) for times much larger than the inverse 
of the characteristic frequency Oo, yet smaller than the 
inverse of the relaxation rate T, is therefore linear in time 
with a characteristic time r + , where 



i , nn 

t+ = — tann — — - . 



(59) 



Note that Eq. ( p9| ) holds for finite temperatures as well 
as in the quantum limit. The phase diffusion time r + is 
inversely proportional to T at high temperatures, 



1 ?if2 
r2fcT' 



(kT > ftfi„), 



(60) 



just as the other characteristic time r_. In the low tem- 
perature or quantum limit, however, it saturates to a 
value 



(61) 



see also Fig. |[ The crossover from high temperature 
behavior to the quantum limit behavior takes place at 
kT ~ hilo. We point out that r + is related to the spec- 
tral densities S^iui) and Sgg(u>) of the internal phases 
(f and 9. This indicates that there is a relation between 
dephasing (at finite temperature) and the reduction of 
the persistent current at zero temperature which is de- 
termined by S vv (lo). 



1 



0.1 



0.01 




2kT/fiOo 



FIG. 5. The characteristic phase diffusion time r+ normal- 
ized by the decay rate F, defined in Eq. j23|), as a function 
of 2kT/hQ,o in a logarithmic scale. For high temperatures, 
r+ is inversely proportional to the temperature T. For low 
temperatures, kT <C HQq, it saturates at a value r + = 



We point out that Eq. ( J55[ ) and Eq. (|56|) do not hold 
for arbitrarily long times. In reality the mean square 
displacements (\c±(t) — c±(0)| 2 ) are bounded since the 
wave function tpu(t) is normalized to 1. The fact that 
(\c-(t) — c_(0)| 2 ), Eq. (p5|), grows without bounds is 
an artifact of the linearization of Eqs. (|l^)-(|l6|) and 
Eq. ©. 

In summary, we find two characteristic phase diffusion 
times t_ and t+, related to the projection of the equi- 
librium state ip(t) onto the ground state and the excited 
state. The time r_ , associated with the projection on 
the ground state, is related to the loss of coherence in the 
global phase \ and diverges as the temperature goes to 
zero. The time t + , on the other hand, characterizes the 
loss of coherence in the internal phases ip and 8, and sat- 
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urates to a finite value in the limit of zero temperature. 
This indicates that even at zero temperature, coupling 
to the external circuit causes excitations of the ring into 
the excited state. In the zero-temperature limit these 
excitations still decay in a finite time. 



V. CONCLUSION 

We have investigated the persistent current in a normal 
metal ring coupled to a resistive external circuit using a 
Langevin equation approach. We have shown that the 
quantum fluctuations in the external circuit suppress the 
persistent current at zero temperature, thus confirming 
an earlier Bethe ansatz result .tH Within the same frame- 
work, but at finite temperature, we have derived two de- 
phasing times t_ and r + , which are the phase diffusion 
times of the projections of the wave function of the ring 
to ground state and the excited state, respectively, of a 
ring which is disconnected from the external circuit. We 
show that t + is related to the spectral densities of the 
"internal" phases tp and 9 of the wave function of the 
ring, which are also responsible for the reduction of the 
persistent current. While r_ diverges in the zero temper- 
ature limit, t + saturates to a finite value. 



APPENDIX A: COULOMB INTERACTIONS 
AND DISPLACEMENT CURRENTS 

In this Appendix, we briefly discuss the electrical cou- 
pling of the ring and the external circuit. The fluctua- 
tions of the charge in the ring and the external circuit are 
not arbitrary but are connected by the flow of a displace- 
ment current between the ring and the external circuit. 
We can express the coupling between the ring and exter- 
nal circuit in terms of the displacement current. Since the 
coupling between these subsystems is crucial, we present 
first a discussion of the displacement currents. 

In our model the charges and potentials are related via 

Qo = C 1 (V -U i ), (Al) 

Qd = C 1 (U d -V )+C i (U d -U a ), (A2) 

Qa = C 2 (U a -Voo)+Ci (U a - U d ) , (A3) 

Qoo = CafVoo - U a ), (A4) 

where we have introduced the parallel capacitance Cj = 
Cl + Cr. In the following, we shall also need the exter- 
nal (serial) capacitance C^ 1 = C-f 1 + C 2 ~ 1 , as well as the 
parallel and the serial total capacitances C = Ci + C e 
and Cq 1 = C,- + C^ 1 . The whole structure, that is, 
the ring together with the external circuit is c harg e neu- 
tral. As a matter of fact, it follows from Eqs. (A1)-(A4) 
that Qo + Qd + Qa + Qoo = 0. The ring and the ex- 
ternal circuit taken separately, however, do not need to 
be neutral. The total charge of the ring Qd + Q a is bal- 
anced by the charge on the external capacitors C\ and Ci , 



Qo + Qoo = -(Qd + Qa), as follows from Eqs. QA1])-(|A4|). 
As the ring and the external circuit do not exchange par- 
ticles, the charges on the ring and in the external circuit 
are conserved, and it is more convenient to consider the 
deviations from a reference state. We denote these devi- 
ations by SQd on the dot and by SQ a on the arm, and 
similarly by SQq and SQoc for the external circuit. They 
obey the relations 5Qd + SQ a — SQo + SQoo = 0. For time 
derivatives, it is of course irrelevant whether we consider 
total charges or deviations. 

The current flowing out from the capacitor C\ is a pure 
displacement current 



h = -Qo = C 1 ^- t (U d -V ), 



(A5) 



whereas the currents flowing through the right and the 
left junction, respectively, are particle currents l v L j R aug- 
mented by displacement currents 



iL = r p L + c L ^- t (Ud-u a ), 
iR = i p R + c R ^- t (u d -u a ). 



(A6) 
(A7) 



The particle currents are related to the charge on the dot 
by /£ + I R = —Qd- At each node we have h + Il + Ir 
I j — (I L + I R ) — 0, and thus I\ + h = 0. These equa- 
tions correspond to the law of current conservation at the 
nodes of an electrical network. 

We can express all electrical quantities of interest in 
terms of the external bias SV = Vo — Voo an d th e charge 
deviation, SQd- We combine Eqs. ( |A2^ ) and (A3) and the 
condition SQd + SQ a — to obtain the potential differ- 
ence inside the ring 



TT TT Ce XW 1 

U d -U a = -gOV+ — . 



(A8) 



The charg e de vi ation SQq on the capacitor C± is found 
from Eqs. (00), (|X|) and (|a|), 



SQ = C SV - ^SQ d , 

where we have used the identity 
C e _ Co 

~c ~ q: 



(A9) 



(A10) 



whence we obtain for the current flowing in the external 
circuit 



h = ^Qd - C SV. 



(All) 



In the circuit containing a current noise source, see Fig. |^, 
the current I\ flowing out of the circuit augmented by the 
current SI coming from the noise source must equal the 
current through the resistor / = 5V/ R, namely 

SV/R = h + 51. (A12) 

Together with Eq. flAll| ), this gives Eq. (0). 
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APPENDIX B: THE CIRCULATING CURRENT 

The equilibrium persistent current I is a quantum and 
statistical average which can be obtained from the deriva- 
tive of the free energy I = —cdF/d<&. Some discussion 
is required, if we are, as in Ref. |l| concerned with current 
fluctuations. Naively, we might want to investigate the 
fluctuations of the persistent current by considering the 
second order derivative of the thermodynamic potential. 
However, such a procedure works only if the observable 
of interest commutes with the Hamiltonian. Moreover, 
as pointed out above, the true and physically relevant 
currents is the total current (particle current plus dis- 
placement current). Whereas the displacement current 
needs not to be considered as long as we are interested 
in average quantities only, this is not true, if we consider 
fluctuations of the current. For the model considered 
here, we can derive expressions for the current operators 
which are particle operators weighted according to the 
distribution of the geometrical capacitances of the sys- 
tem. 

Consider first for a moment the electrically isolated 
loop. In this case the capacit ance s to the exterior circuit 
are C x = C 2 = 0. From Eq. @ we find Q d = d{U d - 
Ua) and thus we can also writ e th e curr ents through the 
left and right junctions, Eqs. (A6) and (A7), in terms of 



the charge on the dot as, II = !'[ + {CL/Ci)dQd/dt and 



Ir 



- TP 



(Cn/Ci)dQd/dt. Using particle conservation, 



L t ± d — dQd/dt = 0, to eliminate the time-derivative 



11+ II 

of the charge we find immediately that the current cir- 
culating in this loop is, I c = II = —Ir with 



Ir, = 



CRlj-Cj R 

Ci 



(Bl) 



where d = Cl + Cr. The circulating current is thus 
in general not determined by the particle currents but 
by an average of these currents weighted according to 
the Coulomb interaction (capacitance ratios). Expres- 
sions of this type are familiar form the dynamic trans- 
port through double barriers, but seem to be novel for 
persistent currents. We emphasize that Eq. (Bl) does 
not mean that the particle currents can now be calcu- 
lated from the non-interacting problem. The dynamic 
particle currents depend on the self-consistent potential 
distribution. For an illustration of this point we refer the 
reader to Ref. |3^ where dynamic current noise spectra for 
double barriers are compared based on calculations using 
the particle currents of the non-interacting problem and 
calculations using the particle currents of the interacting 
problem. 

Next let us consider the system coupled to the ex- 
ternal circuit. In this case C\ > and C2 > 0. We 
can proceed as above. We first express the time deriva- 
tive of potential difference Ud — U a in terms of the time 
derivative of the charge on the dot and the current 



h = dQo/dt = Cid{V Q - U d )/dt (see Eq. (gTj)). Using 
current conservation, I\ = II + Ir and particle conser- 
vation 1 V L + I ft — dQd/dt = we can again eliminate the 
time derivative of the charge Qd and find for the currents 
through the left and right barrier, 



I L = I C + (C L /Ci)h, 



(B2) 



(B3) 



with I c as given in Eq. (Bl). The current I\, which can 



be induced with the help of the external circuit, is di- 
vided onto the two branches of the ring-dot sample also 
according to a capacitance ratio. 

What is observed in a measurement of the magnetiza- 
tion? It is important to note, that the external circuit 
(see Fig. |l|) also forms a current loop and depending on 
the geometry of the circuit also contributes to the fluc- 
tuation of the magnetization. Suppose, the circuit of 
Fig. [I] is a planar structure in the x-y-plane. Let the ring 
with the in-line dot enclose an area A r and let the exter- 
nal circuit (excluding the ring) enclose an area A e . For 
this circuit the magnetic moment can be viewed as being 
generated by a current II enclosing the area A\ and by a 
current I\ enclosing the area A e . The magnetic moment 
is thus m = (l/c)(JiA r . — IiA e ). Of course the topol- 
ogy of the current distribution matters: If we consider 
the external circuit to be above the ring dot structure 
instead of below as shown in Fig. [I], the magnetization 
is m = (l/c)(—lRA r + IiA e ). In these two circuits the 
external circuit gives magnetization contribution into op- 
posite directions. The average magnetization of the two 
circuits is just m = (l/c)I c A e . A third circuit which per- 
mits to investigate the magnetization of the ring and the 
external circuit separately is a structure in which the ring 
lies in the x-y-plane and the external circuit say in the 
s-z-plane. Then the magnetization is m z — (l/c)I c A r 
and m y = (l/c)IiA e . I c appears as the most natural 
generalization of the equilibrium persistent current. In 
Ref. [I] we have investigated the fluctuations of I c . 

All the above expressions are based on linear relations 
between currents, charges and potentials. All the above 
expressions are thus valid also for operators. We now 
present a specific expression for the operator of the cir- 
culating current used in Ref. [j]. Let us consider the parti- 
cle current operators I V L through the left and I V R through 
the right tunnel barrier. In terms of the Hamilton opera- 
tors Hl and Hr responsible for the tunneling across the 
left and the right tunnel barrier, respectively, they read 
l/r = ~i/^\IlL/Ri Qd]- In other words, they are equal 
to the decrease of the charge on the dot per unit of time 
through the right and left junction. For the two-level 
system, in a basis in which the Hamiltonian is real (see 
Eq. (O)) we have 
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e ±t L t R sin(27r$/$ ) a x + [t\ ± t L t R cos(27r$/$ )]g, 

h y/tl + ^±2i L t fl cos(27r$/$ ) 

e ±t^tg sin(27r$/$ ) cr^ - [t| ± tjjg cos(27r$/$ )]g,, 
ft V^l + *1 ± 2t L t R cos(2^$/$ ) 



Using the particle current operators in Eq. (Bl) leads to 
the operator given by Eq. (Q) . With the help of this op- 
erator we can investigate the average persistent current 
and the fluctuations (see Eq. (7) of Ref. |). 
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